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Chandra Measurements of a Complete Sample of X-ray 
Luminous Galaxy Clusters: the Gas Mass Fraction 
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1 INTRODUCTION 



ABSTRACT 

We present Chandra X-ray measurements of the gas mass fraction out to r5oo for a 
complete sample of the 35 most luminous clusters from the Brightest Cluster Sample 
and the Extended Brightest Cluster Sample at redshift z — 0.15 — 0.30. The sample 
includes relaxed and unrelaxed clusters, and the data were analysed independently 
using two pipelines and two different models for the gas density and temperature. 
We measure an average of /gas0"50o) = 0.163 ± 0.032, which is in agreement with the 
cosmic baryon fraction (f2b/^M = 0.167 ± 0.006) at the la level, after adding the 
stellar baryon fraction. Earlier studies reported gas mass fractions significantly lower 
than the cosmic baryon fraction at rsoo, and in some cases higher values that are 
consistent with the cosmic baryon fraction towards the virial radius. In this paper we 
show that the most X-ray luminous clusters in the redshift range z = 0.15 — 0.30 have 
a gas mass fraction that is consistent with the cosmic value at rsoo- 

Sun et all [20091; lUmetsu et all 120091 ; iRasheed et ail |20ld : 



Clusters of galaxies are the largest known bound systems in 
the Universe and are formed from the collapse of primordial 
densi ty fluctuations l|Press fc S chcchter 19741; Iwhite fc Reesl 
Il978h . The intracluster medium (ICM) is in the form of 
a hot ionized gas at ~ 10 s K, and it contains most of 
the baryons in clusters. The remaining baryons are in 
stars and intracluster light, and account for a few per- 
cent of the total mass i|Lin et al.ll2003l ; iGonzalez et~ai1l2007l ; 
iGiodini et al.l 12009). The total cluster baryon fraction is 
therefore a combination of the baryons in stars and ICM, 
i.e., / b = /stars + /gas- Since clusters are extremely large 
and massive, baryons and dark matter originated from ap- 
proximately the same comoving volume, and thus it is be- 
lieved that their ratio should be rep resentative of the Uni- 
verse (e.g., iMetzler fc Evrardl Il994h . A number of stud- 
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evolution of the gas mass fraction over a large redshift range, 
/gas (2), allows cos mological quantities to be constraine d (for 
recent reviews see I Allen et al.ll2008l ; lEttori et al]|2009l and 
references therein). 

Current studies indicate that the cluster baryon frac- 
tion is typically lower than the cos mic baryon fraction as 
measured by the f2b/flM parameter ( Vikhlinin et al.l | 2006| ; 
lAfshordi et al.l2007l ; lArnaud et al.ll2007l ; lGiodini et alj2009t 



Komatsu et al.lfeoill h This observation has r aised questions 
on th e whereabouts of these missing baryons (Rashecd et al. 
|2010| ). A possible explanation is that Active Galactic Nuclei 
(AGN ) feedback may push the ICM towards the cluster out- 
skirts (IMetzler fc Evrardll 19941 ; iTakizawa fc Mineshige] | l998| ; 
iBialek et al.l l200ll; IValdarninil [20031; iMcCarthv et al.l I2007I : 
ICavaliere fc Lapill2008l ; iBode et al.ll2009r i. 

To address this issue we measured the gas mass frac- 
tion for a complete sample of massive clusters at z = 
0.15 — .30 from the Brightest Cluster Sample and its ex- 
tension (|Ebeling et all 1 19981 . 12OO0I ; lDahldl200rj ). The back- 
ground of the surface brightness and temperature profiles 
limits the radius out to which masses can be measured ac- 
curately. Therefore we chose to limit our measurements for 
the entire sample to rsoo, the radius within which the mass 
density is 500 times the critical density of the universe at the 
cluster's redshift. In cases of long exposures and high quality 
data, it is possible to use Chandra to d e tect c lusters beyond 
rsoo- For example, in lBonamente et al.l ()2012h we report the 
detection of A1835, one of the clusters in this sample, out to 
the virial radius. We calculate the total cluster baryon frac- 
tion by adding the baryons from stars and galaxies, and find 
that the baryon content in these high-luminosity clusters is 
consistent with the cosmic ratio fib/^M at rsoo- 

This paper is structured in the following way: Section [2] 
briefly describes the sample of clusters, Section [3] discusses 
the Chandra data reduction, and Section[4]reviews the mod- 
elling of the X-ray data. We present our results in Section [SI 
possible systematic effects in Section [(J a comparison with 
other studies in Section [7] and discussion and conclusions 
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in Section [8] Throughout this paper we assume a cosmology 
based on WMAP7 results with Hp = 70.2 km s " 1 Mpc -1 , 
tt M = 0.27, and fi A = 0.73 (|Komatsu et al.ll201ll ). 



2 CLUSTER SAMPLE 

The sample we chose for t his analysis is composed of 35 clus- 
ters from the lDahld (|2006l ) sample (Tabled), a 90% complete 
sample of clusters with X-ray luminosities in the 0.1 — 2.4 
keV band of Lx.kcV 6 x 10 44 erg s _1 (for a concordance 
ACDM universe w ith h = 0.7) from the Brightest Clus- 
ter Sample f BCS. lEbeling et al.|ll998l ') and the Extended 
BCS (eBCS, lEbeling et al.ll200ol ) in the z = 0.15 - 0.30 
redshift range. These clusters are estimated to have masses 
of Mi8o 5 x 10 14 hT 1 M p), making t hem the most mas- 
sive in the (e)BCS sample (Dahlc 2006). All of the clusters 
from this sample have archival Chandra data available for 
our analysis. The majority of the clusters have been ob- 
served with the ACIS imaging array (ACIS-I), and 6 with 
the ACIS spectroscopic array (ACIS-S). With this sample 
of clusters, we aim to measure the gas mass fraction out 
to rsoo using high S/N X-ray data from Chandra, and then 
calcu late the total baryon fr acti on using the stellar f raction 
from lGonzalez et al.l (| 20071 ) and lGiodini et al.l (|2009l ). 



3 CHANDRA OBSERVATIONS AND DATA 
REDUCTION 

Data reduction was done using the Chandra Interactive 
Analysis of Observations (ciacu) software version 4.2 and 
the Chandra Calibration Database (CALDrQ) version 4.3.1. 
The data reduction and subsequent analysis was performed 
by two separate pipelines: one was developed by D. Landry 
and M. Bonamente, and makes use of a Markov chain 
Mont e Car lo method of analysis (e.g., iBonamente et al] 
|2004 120061 ). the other was developed by P. Giles and 
B. M aughan (additional details are provided in lGiles et alj 
|2012| ). Results from the two analyses, including tempera- 
ture profiles and mass measurements, were found to be in 
statistical agreement, thus providing confidence on the re- 
sults provided in t his paper and in o thers to follow using the 
same sample fe.g.. lGiles et al.ll2012l ). 

As part of the data reduction, corrections were made 
for afterglows, charge transfer inefficiency (CTI), bad pix- 
els and solar flares. Afterglows are caused from cosmic rays 
building up charge on the CCD, bad pixels take into account 
hot pixels and afterglow events, and CTI is due to proton 
damage to the ACIS chips that reduces the energy resolu- 
tion. The variability of solar activity can cause periods of 
high background which need to be filtered out. A common 
way of removing these periods o f solar flares is to follow 
the lightcurve filtering method of iMarkevitch et al.l l|2003l ). 
Before the lightcurve can be created, point sources of high 
and variable emission need to be excluded. The lightcurve 
is then created over the energy range 0.3 — 12.0 keV for 
ACIS-I observations and 2.5 - 6.0 keV for ACIS-S obser- 
vations on a selected region of the Chandra CCD's used 



1 http://cxc.harvard.edu/ciao/ 

2 http://cxc.harvard.edu/caldb/ 
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Figure 1. A sample lightcurve for observation 9371 of 
Zwicky 3146, showing excluded time intervals due to solar flares 
(red boxes). Lightcurves are created over the whole energy range, 
0.3 — 12.0 keV, using the local background region, 450 — 900" in 
this case. The mean count rate is 0.596 s — 1 (dashed line) and the 
filtered exposure time is 36.3 ks. 



as the local background. A sample lightcurve is shown in 
Figure [T] This lightcurve was filtered with an iterative algo- 
rithm (def lare command in CIAO) that removes time inter- 
vals outside the 3a range of the mean. The solid line shows 
the mean count rate for the observation and the red boxes 
show the regions that were filtered out. Most observations 
were taken in VFAINT mode, and in this case we applied 
VFAINT cleaning to both the cluster and blank-sky obser- 
vations. 

Images were created in the 0.7 — 7.0 keV band to mea- 
sure the surface brightness as a function of radius. This en- 
ergy band was chosen since in this energy range the effective 
area is highest and the background rate lowest. The most 
crucial aspect of the analysis of diffuse sources such as galaxy 
clusters is background subtraction. For this purpose, we use 
ACIS blank-sky composite event files and measurements of 
the local backgroun d taken from source-free regions of the 
cluster observation (Ma rkevitch et all 12003). Since all the 
clusters in the sample have redshifts in the 2 = 0.15 — 0.30 
range, the emission does not extend across the entire detec- 
tor, and we have a sufficiently large region to obtain a local 
background (see Figure [2). In nearly every cluster the local 
background was measured from an outer annulus beyond 
»"200, where the surface brightness was approximately con- 
stant, i.e., it has reached the background level. There are a 
few clusters that only have ACIS-S data available. For these 
clusters, the local background was taken from the adjacent 
ACIS-I chips. Point sources and extended substructures were 
detected and removed using wavdetect in CIAO and known 
sources were removed manually by visual inspection. 

The local background of the cluster observation may dif- 
fer from the ACIS blank-sky composites, since these obser- 
vations were done at different times and positions in the sky, 
and the X-ray backgr ound has both spatial depe ndence (e.g., 
Snowden et aLlll997h and is time-variable fe.g.. iTakei et all 



200fet). Although the b ackgr ound flux may vary with time, 



Hickox &; Markevitchl (|2006l ) have shown that the ratio of 
the flux within the 2 — 7 keV and 9.5 — 12 keV bands is 
constant in time. We can therefore rescale the blank-sky 
spectrum by the ratio of the count rates in the 9.5 — 12 keV 
band of cluster and blank-sky observations, and obtain a 
clean subtraction of the background in the spectral region 
of in terest, which is 0.7 - 7.0 keV jrlickox fc Markevitchl 
2006). After subtracting the blank-sky background, resid- 
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Table 1. Cluster Sample 



Cluster 


2 




Ytt 

1 v H 

(10 20 cm -2 ) 


obsID 


Exposure 
fks) 


Dynamical State 


A115 


0.1971 


673.9 


5.36 


3233 


43.6 


Unrelaxed 


A1423 


0.2130 


716.1 


1.81 


538 


35.1 


Unrelaxed 


A1576 


0.2790 


876.1 


1.08 


7938 


15.0 


Unrelaxed 


A1682 


0.2260 


749.5 


1.04 


11725 


19.6 


Unrelaxed 


A1758 


0.2790 


876.1 


1.03 


7710 


7.0 


Unrelaxed 


A1763 


0.2230 


741.9 


0.82 


3591 


17.0 


Unrelaxed 


A1835 


0.2532 


816.3 


2.04 


6880 6881 7370 


167.1 


Relaxed 


A1914 


0.1712 


601.8 


1.06 


3593 542 


16.6 


Unrelaxed 


A2111 


0.2290 


757.1 


1.84 


11726 544 


30.9 


Unrelaxed 


A2204 


0.1520 


545.5 


5.67 


7940 


70.5 


Relaxed 


A2219 


0.2256 


748.5 


1.76 


896 


41.5 


Unrelaxed 


A2261 


0.2240 


744.4 


3.19 


5007 550 


22.5 


Relaxed 


A2390 


0.2329 


766.8 


6.21 


4193 


19.7 


Relaxed 


A2552 


0.3017 


925.9 


4.60 


11730 3288 


22.7 


Unrelaxed 


A2631 


0.2780 


873.8 


3.55 


11728 3248 


25.0 


Unrelaxed 


A267 


0.2310 


762.1 


2.75 


1448 


7.4 


Unrelaxed 


A520 


0.1990 


679.0 


5.65 


528 4215 9424 9426 9430 


368.0 


Unrelaxed 


A586 


0.1710 


601.2 


4.89 


11723 530 


15.8 


Unrelaxed 


A611 


0.2880 


896.1 


4.46 


3194 


15.4 


Relaxed 


A665 


0.1819 


632.1 


4.32 


3586 


14.3 


Unrelaxed 


A68 


0.2546 


819.7 


4.96 


3250 


9.2 


Unrelaxed 


A697 


0.2820 


882.8 


2.93 


4217 


15.4 


Unrelaxed 


A773 


0.2170 


726.5 


1.28 


3588 5006 533 


37.3 


Unrelaxed 


A781 


0.2987 


919.4 


1.65 


534 


9.9 


Unrelaxed 


A963 


0.2060 


697.7 


1.25 


903 


23.0 


Relaxed 


MS 1455+2232 


0.2578 


827.2 


3.18 


4192 543 


79.2 


Relaxed 


RX J0437. 1+0043 


0.2850 


889.5 


5.50 


11729 7900 


41.7 


Relaxed 


RX J0439.0+0715 


0.2300 


759.6 


9.18 


1449 3583 


23.7 


Relaxed 


RX J1720.1+2638 


0.1640 


581.0 


3.36 


3224 4361 


41.7 


Relaxed 


RX J2129.6+0005 


0.2350 


772.1 


3.63 


552 9370 


35.9 


Relaxed 


Z2089 


0.2347 


771.3 


2.86 


10463 7897 


39.4 


Relaxed 


Z3146 


0.2906 


901.8 


2.46 


909 9371 


81.3 


Relaxed 


Z5247 


0.2300 


759.6 


1.61 


11727 539 


19.2 


Unrelaxed 


Z5768 


0.2660 


846.4 


1.49 


7898 


10.4 


Unrelaxed 


Z7215 


0.2897 


899.9 


1.40 


7899 


13.0 


Unrelaxed 




\ mm^ 

2 4 6 8 10 
Right Ascension (J2000) 



Figure 2. The local background region used for ACIST observa- 
tion 9371 of Zwicky 3146- The annulus used for the local back- 
ground is from 450 — 900" (solid lines). The dashed line indicates 
7*500 ar >d is approximately 270". 



uals may still be present in the soft 0.7 — 2 keV band. 
These soft residuals may be due to Galactic and extragalac- 
tic emission, as well as residual sola r flares that were no t 
removed by lightcurve filtering (e.g.. ISnowden et al.lll997h . 
When present, the soft residual spectrum is then fit with a 
power law and a plasma emission model, and this model of 
the soft residuals is taken into account in the cluster spec - 
tra (as done in, e.g.. lBulbul et al.ll20l61 ; lHasler et alj|2012r i. 
Figure [3] shows an example of soft residuals below 2 keV for 
observation 6880 of Abell 1835. Multiple observations were 
reduced individually to apply the correct calibration to each 
dataset. The cluster surface brightness profile is obtained 
from merged images, and the temperature profile from fit- 
ting spectra from different observations simultaneously. 



4 MODELLING AND ANALYSIS OF THE 
CHANDRA DATA 

The total mass of a cluster can be inferred from the density 
and temperature of the X-ray emitting ionized plasma, as- 
suming that the hot gas is in hydrostatic equilibrium with 
the gravitational cluster potential. Thus, we can write the 
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generalisation of the /3-model, 



Energy (keV) 

Figure 3. Spectrum of the local background (beyond ~ 750") 
for observation 6880 of Abell 1835. This shows the soft residuals 
(below 2 keV) fit with an unabsorbed thermal plasma model with 
kT ~ 0.25 keV. 



total mass as 



M tot (r) = 



-fc B T(r) 
firripG 



dlnn c (r) + d In T(r) 
d In r d In r 



(1) 



where fee is the Boltzmann constant, T(r) is the tempera- 
ture profile, fi is the mean molecular weight, m p is the mass 
of a proton, G is the gravitational constant, and n e (r) is 
the number density of electrons. The mass of the hot gas is 
calculated as 



Mga S (r) = 47r/j e w-p / n c (r)r dr 



(2) 



where /i e is the mean molecular weight of the electrons. We 
use high resolution X-ray imaging and spectroscopy from the 
Chandra X-ray Observatory to determine n e (r) and T(r). 
Specifically, the observed X-ray surface brightness is related 
to the electron number density by the following equation 



Sx = 



47r(l + z) 3 



nl{r)K ee {T e ) dl 



(3) 



where z is the redshift, A CC (T C ) is the plasma emissivity 
(counts cm 3 s _1 ) and the integral is along the line of sight 
through the cluster. The linear distance r in these equa- 
tions is given by r = 8 Da, where 8 is the apparent an gular 
size and Da is the angular diameter (|Carroll et al.ll 19921 ). To 
measure the density from the surface brightness, the angular 
diameter must be known. The value of Da can be calculated 
from the following equation (valid for fi^ =0): 



(l + z)H J E{z 



dz, 



where E 2 (z) = Qm(1 + Z ) A + and Ho is the Hubble con- 
stant. The temperature profile is determined from radially- 
averaged X- ray spectra f i t to a n APEC optically-thin emis- 
sion model l|Smith et al.l 1200 ll ) , with Galactic H I column 
density measured fro m the Leiden-Argentine-Bonn survey 
jKalberla et al.|[ 2005). We used the xspec software version 

12.6.0s (|Arnaudll 19961). 

We use the lVikhlinin et al. (2006) model that describes 
the density and temperature of the cluster to fit the Chandra 
data. The three-dimensional gas density is modelled as a 



n p n c 



n * (r/r c )-« 1 , 

(l + ^/r 2 )^-"/ 2 (l + r-y/r]y/-' 

(4) 

n 2 



(l+r 2 /r 2 ) 3 ^ 



and uses a total of ten parameters. The model used for the 
temperature profile is given by the phenomenological func- 
tion 



T(r) 



T (r/r cool ) a ™°' + (T min /T ) 

J 7 — ; s —-. X 

(r/r coo i)«cooi + 1 

(r/n)- a 



(5) 



[l + ( r / rt )6]c/b' 



which has eight parameters and thus enough degrees of free- 
dom to model nearly any smooth temperature distribution. 
The second term in the temperature profile describes the 
region outside of the cool-core as a broken power law with 
a transition region. Therefore, the IVikhlinin et al.l (|2006l ) 
model has a total of 18 parameters. In our analysis, we used 
the following constraints for all clusters: no2 = 0, a = 0, 
7 = 3, and e < 5. For certain clusters with lower S/N, we 
fixed additional parameters. 

To assess biases from the use of a particular parame- 
terization of the th ermodynamic quant ities, we also model 
all clusters with the lBulbul et al.l (|201fj| ) model. This model 
gives analytic functions for temperature, density, and gas 
pressure, assuming a polytropic equation of state and hy- 
dro static equilibrium outs ide of the cluster core. Similar to 
the lVikhlinin et all (|2006t ) model, it accounts for cooling of 
the gas in the core, and the temperature profile is 



T(r) = To 
The function 



1 (l+r/r s f- 2 -1 
(13-2) r/rsil + r/rs)?- 2 

(r/r cool ) a — +^ 



^cool ■ 



"7~cool 



(6) 



(7) 



(r/rcooi^cooi + 1 

has been adopted from I Vikhlinin et al.l (|2006l ) . In Equation 
[3 the quantity £ = Tmm/Tb is a free parameter (0 < £ < 1) 
that measures the amount of central cooling. With this tem- 
perature profile, the gas density profile can be determined 
from the polytropic equation of state, 



n e (r) 



1 (1 + r/n 



N/3-2 



(8) 



,(/3-2) r/r s {l + r/r s y- 2 

The iBulbul et all [j2010h model therefore has a total of 8 
parameters. The comparison of results from the two models 
can be found in Section [6] 

The emissivity for thermal bremsstrahlung emission at 
energy E is 



A e 



-l/2 e -B/fc B T 



where n is the number density and T is the temperature. 
The two X-ray observables, Sx and T, are therefore not 
completely independent. To determine the best-fit param- 
eters of the models, we use a Markov chain Monte Carlo 
(MCMC) method with a Metropolis-Hastin gs selection algo- 
rithm described in lBonamente et al.l l|2004l ). Although both 
models ha ve parameters that can be correlated to one an- 
other (e.g-. lHasler et al.ll2012l ), the MCMC method accounts 
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for these correlations, and model parameters and all deriva- 
tive quantities such as masses and the gas mass fraction can 
be computed accurately. 

In the analysis of the surf ace brightnes s and tem- 
perature profile da ta we follow IBulbul et all (|201Ch and 
lHasler et al.l (|2012h and add a 1% systematic uncertainty 
for each bin of the surface brightness and a 10% systematic 
uncertainty for each bin of the temperature data. In Sec- 
tion [5] we describe additional sources of systematic errors 
that can aff ect our measurements. Selec ted temperature fit s 
for both the lVikhlinin et all (|2006h and lBulbul et all (|2010r i 
model are shown in Figure U (see Appendix |A1 for tempera- 
ture profiles and fits for all clusters). To measure masses at 
rsoo, we need to fit the data beyond this radius to constrain 
the slope of the temperature profile. Out of the entire sample 
of 35 clusters, 28 have temperature and surface brigthness 
measurements out to or beyond rsoo, and only 7 required 
slight extrapolation out to rsoo (see Appendix [X] for addi- 
tional details). 



5 MEASUREMENT OF THE GAS MASS 
FRACTION 

5.1 Results for the Entire Sample 

The gas mass, total mass, and gas mass fraction for each 
cluster are reported in Table [TT] We calculated the median 
and 68.3% confidence interval for the average gas mass frac- 
tion from the combined chains of all clusters, to find 

/ g as(r 2 soo) =0.110 ±0.017 
/ g a S (r 50 o) = 0.163 ± 0.032. 

The average radial profile of the gas mass fraction for all 
clusters is shown in Figure [5] and shows a n increase with 
radius as al so found in previous s tudies, e.g.. lVikhlinin et alj 
(|2006T ) and lRasheed et al.l (|2010l ). 

Since clusters are the largest gravitationally bound 
structures in the Universe, the baryon fraction should be 
representative of the cosmic ratio fib/^M- Current measure- 
ments from WMAP indicate that the cosmic baryonic frac- 
tion is 

n b /n M = o.i67 ±0.006 (10) 

|Komatsu et al.ll201ll ). To accurately compare the gas mass 
fraction with the cosmic baryon fraction, ba ryons in stars 
and g alaxies need to be taken into account. iGiodini et al.l 
(2009) measured the baryon fraction of stars and galaxies 
and determined this stellar mass fraction as a function of 
M 50 o, finding 

/stars, 500 = 0.019 ±0.002 (11) 

for clusters with mass M(rsoo) — 7.1 x 10 14 Mq . This result 
is slig htly higher than the value reported bv lGonzalez et all 
|2007l ). who measu re /stars. 500 — 0.012 for the same mass 
range. Including the lGiodini et all (|2009T ) results in our mea- 
surements of the gas mass fraction, we estimate the baryon 
fraction /b(rsoo) for our sample as 

/b(rsoo) =0.182 ±0.032. (12) 

The difference between the average baryon fraction mea- 
sured at rsoo and the cosmic baryon fraction from WMAP 
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Figure 5. Averag e gas mass fraction profile for the 
Vikhl irTin et alj d2006f ) model for all the clusters in the sample, 
with / gas (r50o) = 0.163 ± 0.032. The red line is the median and 
the hatched region is the 68.3% confidence interval from the com- 
bined Monte Carlo Markov chains for all clusters. The grey en- 
velope is the difference of the cosmic baryon fraction and the 
fraction of baryons in star s and galaxies . Usin g the results from 
iKomatsu et all J201 lT> and IGiodini et al. I d2009h . we use the value 
of U b /£l M - fst ars = 0.148 ±0.006. 

is therefore 

/b(r 5 oo) - Qb/fiM = ±0.015 ± 0.033, 

i.e., our measurements match the known cosmic baryon frac- 
tion at the la level at rsoo, as also shown in Figure [5] The 
measu rement of the baryon fraction using the IBulbul et al.l 
model also agrees with Qb/^M (see Section [6] for dis- 
cussion) . 

5.2 Results for Relaxed and Unrelaxed Clusters 

We used three parameters to determine if a cluster can be 
classified as dynamically relaxed and host a cool-core: cen- 
troid shift for the dynamical state, central cooling time and 
cuspiness for the presence of a cool-core. For clusters with 
high quality data, the central cooling time is th e best method 
for d etermining if it has a cool-core or not (|Hudson et al.l 
2010). We calculated the central cooling time as 

'- = 8.5x10^ (llF^)"(lfe) 1/2 ' (13) 

where n p is the number densi ty of protons and T s is the cen- 
tral temperature of the gas (|Sarazinlll988l ). using the best- 
fitting temperature profile at r = 0.048 rsoo as the central 
temperature. Cuspiness is defined as the logarithmic deriva- 
tive of the density profile, 

a = f"5I* (14) 
a log r 

evaluated at r = 0.04 r 50 o l|Vikhlinin et al.l 120071 '). 

lHudson"et al. (2010) indicate that, for low quality data, this 
parameter is a good indicator for the presence of a cool-core. 
The centroid shift is defined as the standard deviation of 
the di stance between th e peak X-ray emission and the cen- 
troid (|Poole et alj 120061 ) . We measured the centroid shifts 
in annular bins centered on the X-ray peak decreasing by 
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Figure 4. Temperature profiles for selected clusters using the I Vikhlinin et al ] (2006t ) model (blue) and the iBulbul et all l|2010h model 
(red). The solid lines show the best-fit values, and the hatched region is the 68.3% confidence interval. 



0.05 rsoo, as suggested bv lPoole" et all (|2006h . We calculated 
the central cooling times, cuspiness, and centroid shifts for 
all the clusters in the sample, and we classified relaxed clus- 
ters as those which satisfied the following three conditions: 
(i) ( w) < 0.009 rsoo, ( ii) a > 0.65, and (iii) t cool < 6.5 Gyr 
(see iGiles et~ai1l2012l for more details). This classification 
results in 13 relaxed clusters and 22 unrelaxed clusters, as 
shown in Table [T] 

We calculated the gas mass fraction at r2soo and rsoo for 
both subsamples of relaxed clusters and unrelaxed clusters 
(see Table [2j. Our results show that at both radii the gas 
mass fractions for relaxed clusters are in excellent agreement 
with that of unrelaxed clusters, and that both sub-samples 
are statistically consistent with the cosmic baryon fraction 
at rsoo (see Figure |BJ. The agreement of the gas mass frac- 
tion between relaxed and unrelaxed clusters may be an in- 
dication that hydrostatic m ass estimates are r eliable for un- 
relaxed clusters out to rsoo. lGiles et al.l (|2012h discusses fur- 
ther the reliability of hydrostatic mass estimates. 



6 SYSTEMATICS AND PHYSICAL 

PROCESSES AFFECTING THE GAS MASS 
FRACTION 

This section describes certain sources of systematic error 
and physical processes that may affect X-ray measurements 
of the gas mass fraction. 



Table 2. Measurements of /gas using the Vikhlinin et al. (2006 
Model 



Sample 


/gas 




^2500 


^500 


Relaxed 


0.111 ±0.017 


0.155 ±0.026 


Unrelaxed 


0.108 ±0.017 


0.173 ±0.036 


All Clusters 


0.110 ±0.017 


0.163 ±0.032 



6.1 Assumption of hydrostatic equilibrium 

One way in which the gas mass fraction can be affected is 
from the assumption of hydrostatic equilibrium. If there is 
an additional component of pressure that is non-thermal, 
then this must be taken into account when applying the 
equation of hydrostatic equilibrium. The hydrostatic equi- 
librium equation is used to determine the total mass of the 
cluster from X-ray observations, 



M t , 



r 2 dP 



(15) 



where p and P are the density and pressure, respectively, 
of the hot gas. Continued accretion of the gas onto clusters 
along filaments, mergers, and supersonic motions of galaxies 
through the intracluster medium are believed to cause gas 
motio ns which give rise to non-thermal pressure (|Lau et al.l 
2009). Suppose that the pressure of the hot gas consists of a 
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Figure 6. Left: Average gas mass fraction profile using the IVikhlinin et alJ |2006) model for the relaxed clusters, with /gas(r5oo) = 
0.155 ± 0.026. Right: Average gas mass fraction profile for the unrelaxed clusters, with /gasCrsoo) = 0.173 ± 0.036. The red line is the 
median and the hatched region is the 68.3% confidence interval. The blue region is the average / gas profile for all clusters. The grey 
envelope is the diffe r ence of the cosmic ba ryon fraction and the fraction of baryons in stars and galaxies, Sl^/^M — /stars = 0.148 ± 0.006 
jGiodini et a l. 2009; Komatsu et alj|201ll) . 



thermal and a non-thermal component, P - 
Then, the equation for total mass becomes 



-Pth + Pjx 



Af trt (r) = 



r 2 (dPg 



th rf^non-th 

pG \ dr dr 



(16) 



and yields a higher mass compared to the mass from thermal 
pressure only. Separating the mass into components which 
correspond to different pressure terms, the total mass can 
be written as 

Mtot = M th + M non _th. (17) 

lLau et alJ (|2009h shows that the hydrostatic mass underes- 
timates the true mass for several simulated clusters, espe- 
cially at large radii. Gas motions are expected to cause a 
non-thermal pressure in the amount of > 5% — 15% of the 
thermal pressure, and this non-thermal pressure will cause 
an underestimate of mass at large radii of 8% ± 2% for re- 
laxed systems an d 11% ± 6% for unrelaxed systems at rsoo 
(|Lau et al.ll2009l ). In the presence of non-thermal pressure 
the true mass of the cluster is given by Equation 1161 and 
the use of Equation [15] leads to an underestimate of the 
mass , and there f ore an overestimate of the gas mass frac- 
tion. (Giles et al.l (|2012h compared X-ray hydrostatic masses 
with weak-lensing masses for these clusters and found that 
the total mass obtained through X-ray measurements is un- 
derestimated by a factor of 1.21 ± 0.23 and 1.41 ± 0.15 for 
relaxed and unrelaxed clusters, respectively. This compari- 
son would indicate a departure from hydrostatic equilibrium 
that c auses the gas mass fraction to be overestimated. Based 
on the lLau et al.l (|2009h results, we assess a systematic un- 
certainty of +10% on Mtot (f 500 ) for all clusters. This results 
in a possible systematic error of —10% on / ga s(?"5oo)- 

6.2 Uncertainties in the calibration of the 
Chandra data 

Uncertainties in the calibration of the Chandra X-ray data 
can also affect the measurement of the gas mass fraction. 
The efficiency of the ACIS detector has a spatial depen- 
dency, and deviates from being uniform at the ±1% level 



llBulbul et al J2010h . Therefore, we added a ±1% uncertainty 
to each data point in the surface brightness data prior to the 
analysis. 

The temperature profile used in the analysis is also sub- 
ject to various sources of systematic uncertainty. One source 
of uncertainty is from the subtraction of the local back- 
ground. To subtract the background, we first have to rescale 
the blank-sky spectrum to match that of the high-energy 
(9.5 — 12 keV) flux of the cluster observation, as described 
in Section [3] We estimated the effect of the background sub- 
traction using the longest observation of Abell 1835 (ObsID 
6880) which has the highest S/N. To obtain a clean back- 
ground subtraction, we found that the fractional correction 
to the blank-sky spectrum is —0.04 ± 0.01. We measured 
the temperature in an outer region out to ~ rsoo (radii 
240 — 330") using the best-fit correction factor (—0.04) and 
obtained kT = 5.77 keV. We then changed the correction 
to the blank-sky spectrum to —0.03 and —0.05 to over- 
and under-subtract the background by ±1<t. With these val- 
ues, we found the temperature changed to kT = 5.38 keV 
and 6.23 keV, respectively. Therefore, we concluded that 
the uncertainty in temperature due to the background sub- 
traction is approximately ±7.5%. Another source of error 
is caused by contamination on the Optical Blocking Fil- 
ter, w hich is known to affect cluster temp eratures by up 
to 5% l|Bulbul et alj|20ld : lHasler et aTl^Old l . Adding these 
errors in quadrature, we use a ±10% systematic uncertainty 
in fitting the temperature data. This additional error was 
therefore added to each temperature datapoint prior to the 
analysis. 



6.3 Uncertainties due to asphericity of the cluster 
emission 

To measure cluster masses, we assume spherical symmetry 
even though many clusters have a disturbed morphology. 
To estimate the uncertainty due to spherical symmetry we 
considered one of the most disturbed clusters, Abell 520. 
We used observation ID 9426 since this observation had the 
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Figure 7. Image of Abell 520 showing the two sections analysed. 
The regions were s elected by using the temperature map from 
iGovoni et al] \2004) . 
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Figure 8. Temperature profiles for the two sections analysed 
from Figure [7] The blue lines correspond to the eastern section 
and red corresponds to the northern section. The black data is the 
temperature profile obtained from analysing the whole cluster. 



best aimpoint for our task. After performing the Chandra 
data reduction as described in Section [3l this observation 
was left with 96.5 ks of filtered expos ure time. Using the 2 - 
D temperature map of Abell 520 from lGovoni et al.l l|2004) , 
we studied two sections of the cluster (see Figure [7]): the 
northern section (with respect to the azimuthal angle) was 
chosen because it encompasses a gas with temperature ~ 
10 keV, while the eastern section was selected beca use of 
its cooler temperature of < 8 keV jGovoni et aI1l2004l ). The 
two sections are representative of an extreme case of dis- 
turbed dynamical state that causes azimuthal differences in 
temperature and surface brightness. We extracted a tem- 
perature profile and surface brightness profile for the two 
sectors, and measured masses using only the data within 
these azimuthal angles. Figure [H] shows the temperature fits 
to the two sections. For the northern section, we calculated 
/ga a (r 2 50o) = 0.109 ± 0.010 and / gas (r 5 oo) = 0.179 ± 0.015; 
for the eastern section, / ga .s(r25oo) = 0.101 ± 0.007 and 
/gaa(r5oo) = 0.202 ±0.015 (see Table [3] for more informa- 
tion). Using these measurements we estimate that there is a 
±6% and ±8% systematic uncertainty in the gas mass frac- 
tion due to the assumption of spherical symmetry at r2soo 
and rsoo, respectively. 



6.4 Effects of model parameterization 

Another source of error is due to the choice of model 
for the X-ray data. I n addition to the baseline model by 
IVikhlinin et al.l (|2006t ). we r epeated all m e asure ments of the 
gas mass fraction with the iBulbul et al. | (|2010l ) model (see 
Section [4] for more detail about the models). Model param- 
eters for both fits are reported in Tables [SJ [5] and 1 101 (fixed 
parameters are shown without error bars), and the mass 
measurements using both models can be found in Table [TT1 
A plot of the ave rage gas mass fracti on profile for the entire 
sample using the lBulbul et al.l l|2010r i model is given in Fig- 
ure jihowme anjeMellent agreement at all radii with the 
IVikhlinin et all (|2006t ) model results. 

We also compare the average values of ta, M gas , Mtot, 
and /g as for the two models in Table[4j plots comparing each 
measurement can be found in Figure [TOl By using these two 



different models to measure cluster masses we conclude there 
is a 6% ± 2% uncertainty on / gas at r2soo and a 1% ± 2% 
uncertainty on / gas at rsoo due to modelling the data, i.e., 
the two models give the same answer at large radii. 

We also used two independent pipelines to reduce and 
analyse the same Chandra data and measure masses, one de- 
veloped by Landry and Bonamente (LB) and one by Giles 
and Maughan (GM) . The two analyses resulted in the follow- 
ing differences in the measurement of the gas mass fraction: 



/ ga s(r2500) 

/g as (r 5 oo) 



GM-LB 

: GM 
GM - LB 

: GM 



= -2% ± 3% 
= -3% ± 3%. 



This result indicates that our measurements of the gas mass 
fraction using the two pipelines are consistent with one an- 
other. We therefore believe that our measurements of the 
gas mass fraction are robust to the various choices made 

during the analysis. 

We note that a recent study by iMantz fc Alien! (|201lh 
shows how parametric models for the density and tempera- 
ture of the ICM introduce an implicit prior due to the as- 
sumption of hydrostatic equilibrium. If the models are not 
flexible enough or too general, then the deri ved scaling rela- 
tions will be biased towards self-similarity (IMantz fc Alien! 

Eon]). 



6.5 Uncertainties due to clumping of the gas 

Clumping of the gas will also affect the gas mass fraction. 
The main process for X-ray emission in clusters is through 
bremsstrahlung emission, which is proportional to n\ . If the 
gas is clumped, instead of being distributed uniformly, then 
the density of electrons n e will be overestimated by making 
the assumption of a uniform distribution. Observations and 
simulations show that gas clumping may be most evident 
in the outskirts of clusters (r > rsoo, Mathiesen et al.lll999l ; 
iNagai fc Laull201ll ; lEckert et alll2012l ). though even within 
rsoo the g as mass can be overes timated by ~ 10% due to 
clumping (Mathiese n et al]|l999l ). Therefore, our measure- 
ments of /g as at rsoo would remain consistent with the cos- 
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Table 3. Cluster Properties of Sections of Abell 520 



A = 2500 A = 500 

Cluster r A M ga . B M tot / gaB r A M gaB M tat / gaB 

(arcsec) (1O 13 M ) (1O 14 M ) (arcscc) (10 13 Mq) (10 14 Mq) 

A520 166.5 ±1° 3.16 ±§;|1 2.82 ±° f a 0.112 ± <™ 367.1 ±|J 10.71 ±° f> 6.04 ±S;|| 0.177 ±M| 

Hot Slice 168.1 ±}§ ? 3.18 ±g'|| 2.90 ±g || 0.109 ±0 009 375.5 ±\t i 11.60 ±g-|| 6.47 ±g|l 0.179 ±0015 

Cool Slice 165.0 ±f 6 ? B 2.76 ±g|» 2.74 ±«f 0.101 ±^°°% 338.5 ±}« 9.58 ±g : gg 4.74 ±« 3 . 202 ± o:oiS 



T 



Entire Sample (Bulbul 2010 Model) 




[_, 1 , 1 , 1 , 1 , 1 — I 

0.2 0.4 0.6 0.8 1 

r / r 500 

Figure 9. Average gas mass fraction profile for the lBulbul et al] 
(2010) model for all the clusters in the sample with / ga s(f50o) = 
0.161 ± 0.029. The red line is the median and the hatched region 
is the 68.3% confidence interval. The bl u e enve lope is the sample 
average / gas using the IVikhlinin et al.l II2006T ) model. The grey 
envelope is the difference of the cosmic baryon fraction and the 
fractio n of baryons in stars and galaxies, Q h /Qm ~ /stars = 0.148± 
0.006 llGiodini et al.ll2009l ; iKomatsu et al.ll201ll) . 



Table 5. Systematic Uncertainties and Effects on f t 



Source 


Effect on / 
/gas 2500) /g 


^as 

as(r500) 


Chandra Instrument Calibration 






-Surface Brightness 


±1% 




-Temperature 


±10% 




Hydrostatic Equilibrium 
Spherical Symmetry Assumption 
Modelling of X-ray Data 
Clumping of Gas 


-8% 
±6% 
±6% 


-10% 

±8% 
±1% 

-10% 



mological value of Sib / f^M even after accounting for a 10% 
reduction due to possible clumping of the gas. 



ters. V06 measured / gas out to rsoo for 10 relaxed clus- 
ters observed with Chandra spanning a redshift range z = 
0.02 — 0.23, with a weighted average for these 10 clusters of 

/gas, voe = 0.105 ± 0.002. 

The A07 sample consists of 10 relaxed clusters with z ^ 0.15 
observed by XMM-Newton, with a weighted average for the 
10 clusters of 

/gas,A07 = 0.106 ±0.004. 

The sample used in this work are the most luminous clusters 
from the BCS and eBCS in the redshift range 0.15 «S z sC 0.30 
<|Dahlell2006l ). with a weighted average for the 13 relaxed 
clusters of 

,/gas, relaxed 

= 0.150 ± 0.004. 

Notice that the value reported in Table [2] (/ ga s, relaxed = 
0.155 ± 0.026) is the value of the gas mass fraction for the 
average cluster profile, i.e., obtained from the combination 
of all Markov chains. In this section we use the weighted av- 
erage of the 13 measurements for the relaxed clusters, since 
this number can be compared directly with the averages ob- 
tained from the data published in the V06 and A07 papers. 

As an initial comparison, we checked for clusters used in 
our sample that were also analysed by V06 and A07. Clus- 
ters that overlap with our sample were Abell 2390 (anal- 
ysed in V06) and Abell 2204 (analysed in A07). A07 found 
/ g as(r 5 oo) = 0.126 ± 0.013 for Abell 2204, and we mea- 
sure a value of 0.163 ± 0.010. Our temperature profiles do 
no t agree, and compar is on be tween the A07 result shown 
in iPointecouteau et al.l (|2005l ) and our Figure [3] indicate 
that our temperature profile has lower values at large radii. 
Since the temperature of the cluster is the dominant factor 
in calculating hydrostatic masses, the discrepancy in tem- 
perature profiles can be responsible for the different gas 
mass fractions. For Abell 2390, V06 calculates / g as(?~5oo) = 
0.141 ±0.009. Our temperature profile agrees well with that 
of V06, and our measurement / g as(f50o) = 0.131 ± 0.024 is 
also in very good agreement with V06. 



7 COMPARISON WITH PREVIOUS STUDIES 

7.1 Average gas mass fractions from Vikhlinin et 
al. 2006 and Arnaud et al. 2007 

In this section we compare our result to previous studies 
of /gas using hydrostatic mass estimates. We used the re- 
sults f rom the samples of rel axed clusters b y IVikhlinin et all 
l|2006l . hereafter V06) and lArnaud et ail (|2007l , hereafter 
A07), and compare those results with our 13 relaxed clus- 



7.2 Comparison with the Vikhlinin et al. 2006 
results in a matching mass range 

We also compared the masses of our cluster sample with 
those in the V06 and A07 samples. As seen in Figure 111! 
our clusters are generally more massive than those in the 
V06 or A07 samples, yet there is an overlap in the range of 
masses for all three samples. To compare clusters in a similar 
mass range, we grouped clusters in two bins: bin 1 includes 
clusters in the mass range 2 < A/tot < 5 x 1O 14 M0, and 
bin 2 with M to t > 5 x 1O 14 M0. The weighted average of / gas 
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Table 4. Vikhli nin et al, | i2006T ) and lBulbul et alj d2010l) Model Comparison of Sample Average Cluster Properties 



Model 


»"2500 
(arcsec) 


-^gas 
(10 13 M Q ) 


Mtot 
(1O 14 M ) 


/gas 


Vikhlinin et al. (2006) 
Bulbul et al. (2010) 


137.4 ±34.7 
132.4 ±27.7 


2.80 ±1.50 
2.73 ±1.05 


2.53 ± 1.38 
2.36 ±0.96 


0.110 ±0.017 
0.115 ±0.016 


Model 


^500 
(arcsec) 


M gas 

(10 13 M ) 


Mtot 
(1O 14 M ) 


/gas 


Vikhlinin et al. (2006) 
Bulbul et al. (2010) 


292.2 ± 71.1 
301.2 ± 59.4 


8.07 ±3.28 
8.57 ±2.78 


5.05 ± 2.43 
5.58 ±2.26 


0.163 ±0.032 
0.161 ±0.029 



for these bins are reported in Table [6] showing that there is 
a significant difference between our / gas measurements and 
those of V06 and A07, especially for bin 1. Also, contrary 
to V06 and A07, our results do not show an increase in / gas 
with mass between the two mass bins. 

To further investigate this disagreement, we decided to 
analyse all of the V06 clusters that did not require ROSAT 
data for the measurement o f the background, nam ely A133, 
A1413, A383, and A907 i|Vikhlinin et al.l 120061 ). We ex- 
tracted temperature profiles and obtained surface brightness 
profiles for these four clusters, and then measured masses, 
using the same reduction and analysis procedure as for all 
other clusters in this work. Our weighted average for the 
four clusters is 0.114 ± 0.006 (see Table [7}. This is in agree- 
ment with the results from V06, who measure a value of 
0.109 ± 0.003 for these clusters. We therefore conclude that, 
on average, our method of analysis yields statistically con- 
sistent results to those of V06. 

Of the four clusters in this comparison sample, three are 
in agreement at the lcr level. The only cluster we do not find 
statistical agreement with is A1413: we measure / g as(?"50o) = 
0.161 ± 0.011, whereas V06 reports 0.107 ± 0.007. Since the 
temperature we obtain for A1413 at large radii is lower than 
that of V06, we conclude that this is likely the reason for 
the disagreement, since a lower temperature would result 
in a lower total mass and thus a higher gas mass fraction. 
Differences in the Chandra calibration or other aspects of 
the data analysis are likely responsible for the disagreement 
between our results and those of V06 for A1413. We also note 
that A07 also analysed A1413 and measured /gas(?"50o) = 
0.157 ± 0.015, which is in very good agreement with our 
value. 



Vikhlinin et al. (2006) 
Arnaud et al. (2007) 
Relaxed 




10 14 2*10 14 5xl0 14 10 15 

^tot( r 50o) ( M ©) 

Figure 11. f gas - M plot for IVikhlinin et al.l J200fj) , 
I Arnaud et"aH ||2007|) . and the relaxed clusters from this work. 
The mass ranges show overlap between the samples. 

analysis of clusters in the same mass range indicates that the 
selection of clusters based on X-ray luminosity - as in the 
case of our sample - may result in the preferential selection 
of the high- / gas tail of the cluster / gas distribution for a given 
mass. This is not surprising, since both / gas and Lx depend 
on the gas mass content of the cluster. This concl usion is also 
suppo rted by the scaling relations measured bv lGiles et al.l 
(2012), in which we use the same sample of high-Lx clusters 
to find that there is an offset with respect to earlier studies 
that can be explained by a higher gas mass for a fixed total 
mass. 



7.3 A possible luminosity— selection bias for the 
gas mass fraction 

The sample of 35 clusters used in this work was selected as 
the most X-ray luminous in the 0.15 — 0.30 redshift range. 
We calculated the X-ray luminosity in the 0.6 — 9. keV band 
using spectra within the (0.15 — 1) rsoo region (iGiles et alj 
l2012h . We compare the Lx for the five relaxed clusters in 
bin 1 of our sample, and the three clusters (A133, A383, 
and A907) from the V06 sample in the same mass range. 
A133, A383, and A907 were found to have luminosity of 
respectively 5.5 x 10 43 , 2.8 x 10 44 , and 3.8 x 10 44 erg s" 1 , for 
an average of 2.4 x 10 44 erg s _1 , a factor of approximately 
three times lower than the values for the five clusters in 
the same mass range present in our sample (Table [6}. This 



8 DISCUSSION AND CONCLUSIONS 

In this paper we studied the gas mass fraction for the 35 
most lumin ous clusters in the BCS/eBCS at redshift z — 
0.15 — 0.30 (|Dahlei l2006). In accord with earlier studies, we 
find that the gas mass fraction increases with radius, and 
thus the value of / gas depends on the radius used. 

Low gas mass fractions (e.g., ,f ga s ^ 0.1 ) have been ob- 
served in groups of galaxies dSun et al.ll2009f ). likely because 
the gas is not fully bound to the group due to the shallow 
gravitational potential. Recent studies report that the clus- 
ter baryon fraction falls short of the cosmic baryon fraction 
at rs o o even in more massive clusters (e.g., IVikUinin et aL 
2006 ; I Arnaud et all 120071 ; lEttori et al.ll2009l : lRasIieed et all 
20101 ). In our sample of the most X-ray luminous clusters at 
redshift z = 0.15— 0.30, we find a significantly higher average 
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Figure 10. Comparison between measurements using the IVikhlinin et al ] (|2006h (Vikh) and iBulbul et all j2010h (Poly) models. The 
black line is y = x, unrelaxed clusters are shown in red, and the blue data points correspond to relaxed clusters. 



value of the gas mass fraction at rsoo, /gas = 0.163 ± 0.032. 
To find t he total baryon fract ion we used the stellar frac - 
tion from iGiodini et al.l (|2009l) and iGonzalez et al.l l|2007l) . 
Comparing the total baryon fraction with the cosmic baryon 
fraction (Qh/^M = 0.167 ±0.006) we find that the two mea- 
surements agree at the la level, i.e., the gas mass fraction 
at rsoo is in fact consistent with the cosmological value of 



fib/^M, and there are no missing baryons within rsoo in 
the most luminous and massive clusters. A recent study by 
iMiller et al.l {2012) also find a gas mass fraction in agree- 
ment with the cosmic baryon fraction using data from Chan- 
dra, XMM-Newton, and Suzaku. 

One question that remains open is: what happens to 
/ gaa beyond rsoo? As seen in Figure [5] the gas mass fraction 
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Table 6. Cluster Properties for Mass Bin 1 ( 2 < M tot < 5 X 10 14 M ) and Mass Bin 2 (M tot > 5 X 10 14 M ) 



Sample 




/gai 

Bin 1 


i(r 5 oo) 

Bin 2 


Number of Clusters Lx 

Bin 1 Bin 2 Bin 1 (10 44 erg s" 1 ) 


Vikhlinin et al. (2006) 


0.109 ± 0.004 


0.116 ±0.003 


3 


5 


2.4 


Arnaud et al. f 20071 
Relaxed Clusters 


0.111 ± 0.006 
0.158 ± 0.009 


0.131 ±0.009 
0.147 ±0.005 


3 
5 


3 
8 


7.1 


Table 7. Comparison of Clusters with Vikhlinin et al. (2006) 


Cluster 


^500 

(kpc) 


Mtot(r-2500) 

(1O 14 M ) 


/gas(r2500) 


M to t(r-5oo) 
(1O 14 M ) 


/gas(?"50o) 


Weighted Avg. 
/gas(r50o) 








Vikhlinin et al. (2006) 






A133 
A1413 
A383 
A907 


1007 ± 41 
1299 ± 43 
944 ± 32 
1096 ± 30 


1.13 ±0.07 
3.01 ± 0.18 
1.64 ±0.14 
2.21 ± 0.14 


0.067 ± 0.002 
0.094 ± 0.003 
0.092 ± 0.005 
0.091 ±0.003 


3.17 ±0.38 
7.57 ±0.76 
3.06 ±0.31 
4.56 ±0.37 


0.083 ± 0.006 
0.107 ±0.007 
0.124 ±0.007 
0.124 ±0.006 


0.109 ±0.003 








Landry et al. (2012) 






A133 
A1413 
A383 
A907 


1027 ± 48 
1160 ±40 
886 ± 98 
1142 ± 53 


1.11 ± 0.07 
3.55 ± 0.45 
1.40 ±0.46 
2.19 ± 0.12 


0.066 ± 0.002 
0.097 ±0.007 
0.120 ±0.022 
0.104 ±0.003 


3.26 ± 0.46 
5.06 ± 0.52 
2.36 ±0.79 
4.91 ±0.67 


0.080 ± 0.008 
0.161 ±0.011 
0.169 ±0.045 
0.132 ±0.013 


0.114 ±0.006 



increases with radius. A naive extrapolation beyond rsoo will 
give an even higher value for the gas mass fraction. The total 
baryon fraction in this case will be greater than the cosmo- 
logical value, Qh/^M- High values of the gas mass fractions 
towards the virial radius have been reported in recent studies 
based on Suzaku data (e.g.. [Simionescu et al.lfeoill ). How- 
ever, there are two important reasons why simple extrapo- 
lations to large radii of the measured radial trend of / gas (r) 
may not be valid. First, the gas beyond rsoo is not likely to 
be in hydrostatic equilibrium. Non-thermal pressure support 
will become more significant at radii > rsoo and will im- 
pact the measu rement of the total mass (|Nagai et al.ll2007l ; 
lLau et al.|[2009h . An increase of non-thermal pressure sup- 
port would cause an underestimate of total mass, and thus 
an overestimate of /gas - Therefore, the true gas mass fraction 
would be lower than our measured value in the presence of 
non-thermal pressure. The second reason is that clumping of 
the gas may also become more important at large radii. Ob- 
servations and simulations show that the gas clumping factor 
increases beyon d rsoo and will cons i derably bias X-ray mass 
measurements (iNagai fc Laul l201ll ; ISimionescu et al] l201ll ; 
lEckert et al.ll2012T ). As a result of clumping, the true gas 
mass fraction at large radii would be lower than the value 
measured assuming a smooth distribution of the gas. 

We also consider the effect that our sample selection 
criteria could have on the measurement of / ga s- The sample 
presented in this paper comprises clusters selected because 
of their high X-ray flux and luminosity, i.e., our sample pref- 
erentially selects clusters with a high Lx for a given mass. 
Since Lx depends on the gas mass, a luminosity-selected 
sample may have a larger fraction of clusters with high / gas 
than samples selected with other criteria. The extent of any 



such bias would depend on the correlation between the scat- 
ter in /gas and Lx at fixed mass, and the magnitude of the 
intrinsic scatter in / gas - Full account of this effect requires 
modelling of the selection functio n and population scatter of 
the sample, which is presented in lGiles et al.l (|20T2t h For the 
current study, we estimate the intrinsic scatter of / gas in our 
sample as 11%±4%. This may be underestimated if clusters 
from the low / gas tail of the distribution are excluded by the 
luminosity selection of the sample, ff this bias exists, then 
it could pose problems for using the gas mass fraction for 
cosmology, fn fact, clusters used for cosmological applica- 
tions are generally selected on the basis of their high lumi- 
nosity, allowing their study at high redshift. Such clusters 
would be biased towards high / gas , potentially distorting the 
cosmological constraints derived from the / gas (z) tests. The 
implication is that complete (or at least statistically rep- 
resentative) samples should be used for the f sas (z) tests, so 
that the selection bias can be corrected. While these samples 
would necessarily contain morphologically disturbed clus- 
ters, our results show that / gas measurements are not signif- 
icantly affected b y cluster morphology, at least out to rsoo; 
iGiles et al.l (|2012l ) also find an agreement between scaling re- 
lations for relaxed and unrelaxed hydrostatic masses. This 
suggests that complete samples of clusters could be used for 
/ g as(z) tests. A complementary approach would be to study 
the baryon fraction of clusters selected independently of the 
ICM (e.g., red sequence or weak lensing selected clusters). 
This would give a measurement of the range of baryon frac- 
tions that is free from selection biases associated with their 
X-ray emission. 

We conclude that the large value of / gas at rsoo mea- 
sured for this sample is representative of the high-/ gas tail 
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of the cluster population, and that the most massive and 
X-ray luminous clusters in this redshift range have the cos- 
mological ratio of baryons to dark matter even at rsoo ■ 
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Table 8. IVikhlinin et al.l i2006T ) Model Density Parameters 



Cluster 


™e0 


r c 




r s 


e 


™e02 


7 


a 




(10~ 2 cm -3 ) 


(arcscc) 




(arcscc) 











A115 

A1423 

A1576 

A1682 

A1758 

A1763 

A1835 

A1914 

A2111 

A2204 

A2219 

A2261 

A2390 

A2552 

A2631 

A267 

A520 

A586 

A611 

A665 

A68 

A697 

A773 

A781 

A963 

MS 1455.0+2232 

RX J0437. 1+0043 

RX J0439. 0+0715 

RX J1720. 1+2638 

RX J2129. 6+0005 

Z2089 

Z3146 

Z5247 

Z5768 

Z7215 



5.86 ±g 
2.58 ±2 
1.16 ±8 
0.64 ±8 
0.37 ± 
0.81 
9.84 ±Sj 

i 47 +o:os 

0.60 ±g : g| 
30.10 ±};f| 
0.90 ±8 
3.10 ±2 



0.01 

8:o°J 



02 

Hi 

Hi 
W 
..48 
8:V5 



7.64 ± 
1.74 ±[ 
0.57 ± 
0.93 ± 
0.37 ± 

I. 63 ± 
2.14 ±j 

i.oo ±: 

0.82 ± 
0.90 ±5 
0.87 ±g 
0.27 ±2 
2.57 ±2 
7.30 ±2 
4.07 ±2 
3.06 ±8 
7.37 ±2 
6.43 ±g 

II. 21 ±f-|| 

10.49 ±*;i 

0.19 ±° : °i 
0.63 ±»;» 

62 ±g:SS 



4-89 ±g:|jj 
9.85 ± ; s3 
23.33 ±|ii 

2o.i9 ±|:oI 

81.45 ±fH 
oe 77 _|_2.74 

6.69 ±8:ji 

54.80 ±|;J| 
42.10 ± 



3.71 ± 
59.85 ± 
13.44 ± 
3.86 ± 



2. 10 
0.10 

f !iMi 
:B§ 

1 . 10 
0.77 



9.42 ±|;|j 
51.52 ±8;°j 
35.98 ±l;ei 
127.30 ±l l% 
31.30 ±°;|| 

16.59 ±} : }| 
23.04 ±\-%l 
43.14 i?^ 
43.42 ±6;0| 
45.61 ±} |f 
111.40 ±ig 
12.66 ±\-}l 

n.47 ±g;»f 

6.48 ±g 
8.69 ±|-j 

9-38 ±°oh 
fi 79 ,0.8? 
°- ' J ±062 
5.46 ±H| 
6.89 ±° g 

84.60 ±6;|| 
15.49 ± 



I 0.91 
Hi. 90 

L 1. 10 



0.42 ± 
0.44 ± 

0.53 ±; 

0.33 ±g 
0.7 
0.48 ±2 
0.50 ±| 
0.73 ±2 
0.54 ±2 
0.46 ±[J 
0.63 ±g 
0.44 ± 
0.36 ± 
0.29 ± 
0.55 ± 
0.64 ± 
0.85 ± 
0.61 ± 
0.56 ± 
0.41 ± 
0.62 ± 
0.58 ± 
0.62 ± 
0.9 
0.43 ± 
0.45 ± 
0.38 ± 
0.38 ± 
0.48 ± 
0.44 ± 
0.54 ± 



!i:!IJ 



!!:?o 

,.o ! i 

11:111 

0.01 
0.01 

ill 
:§ 
.Jl 
!!:!!i 
!!:!!2 

!i:Hi 



280.80 ±; 



_110.30 
c 110.90 
100.0 
190.80 ±ilH 
100.0 



500.0 
10.0 
133.00 ±g;gg 

500.0 
70.00 ±1111 
145.10 ±J?T2 



_0.2(i 
77°79 



51.41 ± 



5 S3 
3.74 
10 
7.84 



51.24 ±| 

178.90 ±40.74 

100.0 

900.0 

10.0 

426.30 ±\$M 

472.20 ±|1'77 

233.40 ±}$ : §° 

279 30 +159.40 
^cy.JU =t 162 .50 

100.0 

500.0 

176.10 ±JSJ2 

4.96 ±g;|§ 

48.88 ±f-|| 

W8.40 ±|J;1; 

97.82 iijill 

76.7i ±fl;fj 

99.43 ±ll 2 5 7 l° 
96.05 ±l% 7 e 9 

100.0 

100.0 

500.0 



0.0 
3.0 
1 IS +0-07 

0.0 

0.0 
i 12 +0.06 
J-.J-s ^o.06 

0.0 

0. 98 ±°;}| 
2.67 ±g;l| 

i-57 ±8:1? 

2.0 
().() 
().() 
0.0 
5.0 
5.0 
2.98 ±HS 

1 04 +0.83 

1. U4 ± 72 

0.0 
0.0 

3.05 ±g;l| 

1.0 
i « +0.16 
±0.15 



i i o +0.15 
J-- 1 ^ ±0.12 

1.0 
i = fi +0.18 
l.DD ±0.17 

0.0 
().() 
().() 



0.0 


3 





().() 


0.0 


3 





().() 


0.0 


3 





().() 


0.0 


3 





().() 


0.0 


3 





().() 


0.0 


3 





0.0 


0.0 


3 





().() 


0.0 


3 





().() 


0.0 


3 





().() 


0.0 


3 





().() 


0.0 


3 





().() 


0.0 


3 





().() 


0.0 


3 





().() 


0.0 


3 





().() 


0.0 


3 





().() 


0.0 


3 





().() 


0.0 


3 





().() 


0.0 


3 





().() 


0.0 


3 





0.0 


0.0 


3 





().() 


0.0 


3 





0.0 


0.0 


3 





().() 


0.0 


3 





().() 


0.0 


3 





().() 


0.0 


3 





().() 


0.0 


3 





().() 


0.0 


3 





().() 


0.0 


3 





().() 


0.0 


3 





().() 


0.0 


3 





().() 


0.0 


3 





().() 


0.0 


3 





().() 


0.0 


3 





().() 


0.0 


3 





().() 


0.0 


3 





().() 



© 2012 RAS, MNRAS OOO.ITH201 



The Gas Fraction of Luminous Clusters 15 



Table 9. Vikhli nTn et al.l ||2006| ') Model Temperature Parameters 



Cluster 


T 
(keV) 


(keV) 




''"cool 
(arcsec) 


a cool 




n 

(arcsec) 


< 


/ 


b 




c 


X?ot(d-o.f.) 


A115 


21.94 ±l° 7 l 


2.08 ±g' 
4.98 ±!' 


36 
06 


71 81 
222.50 ±78^68 


1.0 




307.60 ±72:12 








2.0 




2.0 


92.50 (128) 


A1423 


7.47 ±J;|| 


50.0 


2.0 




183.80 ±66;q6 





.0 


2.0 




1.0 


56.65 (113) 


A1576 


r\ rtn 1 2 , 48 

9.20 ±i. B l 


9.20 


50.0 


1.0 




/ii 1 rn i_373.90 
411.50 ±260.90 





.0 


1.0 




1.0 


16.80 (74) 


A1682 


8.91 ±f|| 


2.0 




10 a _i_8.94 
16.43 ±| g4 


2.0 




101 T/i 1 71 64 
121.70 ±40^0 





.0 


2.0 




1.0 


52.65 (89) 


A1758 


10.40 ±l;g 


10.40 




10.0 


0.0 




100.0 





.0 


1.0 




().() 


104.51 (59) 


A1763 


897 ±5:03 


8.97 




50.0 


1.0 




1200.0 





.0 


1.0 




1.0 


97.75 (126) 


A1835 


21.98 ±§;|j| 


3.0 




51.01 


1.0 




307.10 ±oQ3g 

487.00 ±227'6Q 





.0 


1.14 ±3 


52 
34 


2.0 


196.64 (83) 


A1914 


i i 00 1 1 91 

11.38 ±{ ;^ 


11.38 




50.0 


1.0 







.0 


1.0 


1.0 


66.00 (109) 


A2111 


/i o _i_ 1 - 86 
8.43 ±1.40 


6.0 




50.0 


1.0 




a an on 1 300 . 20 
469.30 ±i S 4. 8 o 





.0 


2.0 




2.0 


43. 9o (y3J 


A2204 


9.94 ±°;f| 
14.16 ±{;|i 


0.09 ±0 


10 
06 


•7 70 _i_0.69 
'■ ' s ± 0.62 


2.0 




~ Q ~ ~ n 1 320.60 
636.70 ±i 8 7.oo 





.0 


1.94 ±q 


57 
48 


2.0 


153.68 (136) 


A2219 


14.16 




50.0 


1.0 




one t?r, 1 230 50 

305.60 ±120 so 





.0 


1.0 




1.0 


21.66 (104) 


A2261 


11.66 ±i; 3 l 


4.0 




30.0 


0.5 




299.80 ±3-010 








2.0 




2.0 


73.71 (105) 


A2390 


30.11 ± 4 |S 


0.0 




60.0 


0.61 ±q 


12 
13 


1000.0 





.0 


1.0 




4.00 ±111 


87.44 (152) 


A 2 5 5 2 


8.77 ±}1| 

9.36 ±1 24 
k no 1 0.91 
O.UH ±q. S6 


8. 77 




10.0 


0.0 


300.0 





.0 


2.0 




2.0 


61.66 I ( 4) 


A 2 631 


9.36 




10.0 


0.0 




316.80 ±?f|;f° 

c ri 1 3 8 5 . l ) 

385.30 ±143.80 





.0 


1.0 




1.0 


29.95 (74) 


A 267 


5.08 




10.0 


0.0 







.0 


2.0 




2.0 


122.55 (152) 


A 52(1 


<-?a 1 0.55 

8 - 76 ±0.58 
oq _i_ . 76 
S.S3 ±0.74 

9 - 86 ±i:57 


8.76 




3.0 


1.0 




*> ag er\ 1 45 . 52 

346.50 ±37 





.0 


2.0 




2.0 


43.71 (128) 


A58G 


5.0 




60.0 


0.2 




330.30 ±64; 6 | 





.0 


2.0 




2.0 


49.60 (103) 


A611 


5.0 




20.0 


1.0 




300.0 





.0 


2.0 




2.0 


26.61 (47) 


A665 


, | 7 , _2.Sb 
11.(1 ±2.02 


6.41 ±1' 


17 
27 


100.0 


1.0 




eln , n 1 206 40 

510.40 ±155:30 





.0 


2.0 




2.0 


180.08 (179) 


A68 


10.08 ±|;|| 


10.08 




50.0 


1.0 




f>n 1 450 00 

269.60 ±i6i go 








1.0 




1.0 


23.82 (78) 


A697 


13.87 ±| jj 


6.0 




00 00 _i_29. 12 
JO. JO =1=22.40 


2.0 




407.10 ±lg|:?g 





.0 


1.0 




1.0 


29.05 (117) 


A773 


7.73 ±g;f| 
6.53 i}}! 


7.73 




10.0 


0.0 




1000.0 





.0 


1.0 




1.0 


127.68 (94) 


A781 


6.53 




50.0 


0.0 




1000.0 





.0 


1.0 




1.0 


34.04 (75) 


A963 


10.65 ±fg 


3.0 




20.0 


1.0 




98.05 ±1§:1? 





.0 


2.0 




1.0 


76.59 (77) 


MS 1455.0+2232 


4 qi + 0.38 
4.ai ±0.42 


3.87 ±S; 


66 

i!> 

38 


20.0 


2.0 




2000.0 





.0 


2.0 




2.0 


150.82 (102) 


RX J0437. 1+0043 


10.38 ±=;°| 
11.54 ±fH 


3.28 ±{ 


30.0 


1.0 




278.90 ±}SJH 





.0 


2.0 




2.0 


49.92 (57) 


R.X J0439. 0+0715 


3.0 




50.0 


1.0 




190.50 ±gg : *° 

363.80 ±111™ 





.0 


2.0 




2.0 


14.34 (90) 


RX ,71720.1+2638 


13.94 ±l;g 


4.37 ±° 


41 
1 1 


73.39 ±}| : p 
32.80 ±}J;f| 


2.0 







.0 


1.0 




1.0 


163.51 (101) 


RX .72129.6+0005 


9.43 ±111 


4.19 ±g; 


71 

82 


2.0 




141 yn _i_70.04 
141. ±41 80 





.0 


2.0 




1.0 


44.31 (74) 


Z2089 


4.86 ±8:ll 


3.5 




50.0 


0.70 ±[] 


63 
49 


942.70 ±68|'ip 





.0 


2.0 




2.0 


6.93 (43) 


Z3146 


10.39 ±f;i° 


2.0 




18.84 ±f°it 5 


1.0 




265.40 ±p 5 4 f 
241.50 ±?i : gg 





.0 


2.0 




2.0 


98.22 (90) 


Z5247 


7.15 ±H\ 
4.55 ±11? 


7.15 




50.0 


0.0 







.0 


1.0 




1.0 


43.82 (82) 


Z5768 


4.55 




50.0 


0.0 




500.0 





.0 


1.0 




1.0 


29.35 (59) 


Z7215 


8.68 ±i\l 


2.0 




17-1-1 1 1 2 . 2 3 
l/.ll =t 10 23 


2.0 




300.0 





.0 


2.0 




2.0 


18.15 (59) 
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Table lO. lBulbul et alj [|2010h Model Parameters 



Cluster 


™eO 




P + l 


To 


r cool 


t 




X? ot (d.o.f.) 




(10~ 2 cm- 3 ) 


(arcscc) 




(koV) 


(arcscc) 









AllS 

A1423 

A1576 

A1682 

A1758 

A1763 

A1835 

A1914 

A2111 

A2204 

A2219 

A2261 

A2390 

A2552 

A2631 

A267 

A520 

A586 

A611 

A665 

A68 

A697 

A773 

A781 

A963 

MS 1455.0+2232 

RX J0437. 1 + 0043 

RX J0439. 0+0715 

RX J1720. 1 + 2638 

RX .12129. 6+0005 

Z2089 

Z3146 

Z5247 

Z5768 

Z7215 



0.80 ±0-16 

,+ 8iil 
;: as 

Si lid 

iiiiH 

: §? 

o':i 7 , 

!i ,v; 

" 83 

m 

8 7! 

m 

8:48 
8:88 
8:8? 

8 94 
8i8 
8 88 
8 V/ 
ii 88 
!i m 
. Ii 
8:58 

Hi 
8 89 
" 38 
i: *A 
8 89 

V 84 

i:§e 

"3f 
.84 
8:H 



1.81 ± 

0.41 ± 

0.51 ± 

1.17 ± 

4.52 ± 

2.27 ± 
0.90 ± 
4.15 ± 
1.38 ± 
2.84 ± 
1.09 ± 
2.24 ± 
0.76 ± 
1.42 ± 
0.53 ± 
0.92 ± 

1.28 ± 
0.38 ± 
1.14 ± 
0.99 ± 
1.33 ± 
0.36 ± 
1.54 ± 
3.27 ± 
2.61 ± 
0.72 ± 
3.78 ± 
1.99 ± 
5.13 ± 
5.22 ± 
0.24 ± 
0.49 ± 
0.59 ± 



•*.s.ui 
20.77 ±f 
272.60 ±}]"» 
347.10 

41.05 ±f f[ 
14.63 ±KJ 
77.65 ±l\A 
42.98 ±|;« 
22.56 ±\ Ij 

8236 ±igiJs 

18.29 ±gZS 
f),H 2 



155.40 ± 
17.51 ± 
404.00 
57.45 ±_ 
496.00 ± 



± il 1 si° 



6835 HUt 

37.50 ±?r|J 
127.20 ±||;g 
225.90 ±oS q2 



56.17 ± 
1442.00 ± 



60 

§*bo 
"1:99 



7.45 ± 
17.57 ± 
288.90 ± 
18.30 ± 



25.84 ± 
6.20 ± 



.Mr 

,i6 2 §8 

:im 



33 

u 



211 



6.72 ± 
2.08 ± 
4.91 ±_ 

4.0 

10.0 

10.0 
8.99 ±i'8f 



34.45 
6.60 
4.93 
5.19 
3.60 



+1:11 

±?1 3 87 

+l 6 o8 7 



3.40 ±J 
11.08 ± 



2.58 ± 
4.56 ± 
3.42 ± 
3.52 ± 



;f. 1:8 



±: 

± 

± 
± 



"18 

1°0!J 7 

m 1 



10.2 

7.11 
33.81 

3.78 

3.19 
48.98 

2.45 : 

14.53 ±J™> 

4.39 

A -4-1 -Od 

± i,-«, 
13.31 ±§;jf 

10.0 

5.0 

9.07 ±S"5g 



2.27 ± 

3.28 ± 
2.53 ± 
3.95 ±[ 

3.0 
2.27 ±' 

2.29 ± 
2.56 ± 
2.08 ± 
2.39 ± 
2.81 ± 
2.48 ± 
3.95 ± 
2.16 ± 
5.39 ±\ 
3.35 ±§ 

3.0 
3.69 ±8 
2.76 ±jj 
3.13 ±8 

5.35 ±J 

2.36 ±g 



i-iiW 

"69 



8:84 

"% 
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Table 11. Cluster Masses 



(arcscc) 



'"gas 

(10 13 M Q ) 



2500 
Mtot 
(10 14 M ) 



(arcscc) 



A : 

A/gas 

(1O 13 M ) 



500 

-Mtot 
(1O 14 M ) 



Vikh 
Poly 



134.5 ± 
139.7 ± 



13.7 
12.7 
4.0 
4.2 



1.57 
1.68 



I 0.28 
±0.25 

i 0.09 
±0.09 



103 + o: "" r ' 
U.1U3 ±q 004 
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Poly 
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1.48 ±2;1| 
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267.3 ±}2 
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Vikh 
Poly 
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Poly 
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127.4 ±1?;1 

122.0 ±l;| 



140.0 ±S 
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2.23 ± 
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0.28 
1} 29 
0.22 
0.21 



i fifi i 0.26 
±0.27 
1 B4 +0.16 
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3.09 
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4.00 
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3.17 ± 
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2 22 

4 74 1 0.60 
*• '* ±0 61 

300 I 0.20 
- ±0.20 



0.107 ±2 
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101 +0.006 
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12.0 
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0.159 
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Vikh 
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4.9 



4.71 ± 
4.82 ± 



0.21 
21 
(1 25 
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12 1 
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Poly 



115.9 ±S 
109.5 ±1 
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105 +0 016 
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15.7 
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269.3 ±11° 



412.7 ±i 
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5.66 ±g;l 4 
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a qi _i_0.83 
± Q-S1 
7.92 ±g;|g 



3.67 ± 
3.76 ± 
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0.57 
0.48 
0.47 



e.oo ±i;Ig 



6.46 ±f 
5.57 ±? 



0.155 
0.151 



0.118 
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0.121 
0.145 



0.129 
0.143 



_ 0.0 4 9 
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1=0.023 
0.017 

8:8ii 

0.009 
0.037 

8ii 
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Vikh 
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4.32 



0.127 ±ii;l-;j§ 

133 +"- 1 ' 11 
u.J-oo ± n nno 



309.9 ±1 
326.8 ±* 



0:2 15-48 ±i;g| 8.68 ±^;i5 
0; 4 16.60 ±}; 2 1 10.18 ±l°% 



0.178 
0.163 



A7S1 



Vikh 
Poly 



Vikh 
Poly 



146.1 ±9. 



69.7 ±11% 



, 71 +0.20 
z - ' 1 ±0.19 
2.80 ±°;°f 



0.71 
0.68 



+ 0.54 
±0.44 
+ 0.23 
±0.19 



2.38 ±; 



0.57 ±°,ll 
0.53 ±8:?I 



0.119 ±8:888 
ens ±i °°i 
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319.8 ±\ft 
334.0 ±g g 



249.0 ±111 



o.ou ±0.98 

9.95 ±i;ij 



4.98 ±2 
5.67 ±° 



±1.40 

7.37 ±l;f| 
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0.161 
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0.135 



Vikh 
Poly 



153.8 ±|;J 

149.9 ±|;§ 



2.66 

2.55 



4_0.20 
±0.28 
1 0.42 
±0.42 



0.110 ±° °ll 

7 



2.25 ±8:53 0.114 ± 



308.2 ±24 ? 
309.5 ±10 



6.57 ±8 : 1| 3.96 ±g;g| 



0.166 
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Table ITT1 Cluster Masses-Continued 



A = 2500 A = 500 

>"A *fga,s M to t /gas >"A *fgas M to% 

(arcsoc) (10 13 Mq) (10 14 Mq) (arcscc) (10 13 Mq) (10 14 Mq) 



Vikh 
Poly 



Vikh 
Poly 



104.9 ± 
105.3 ± 



3.7 
1.2 

2.5 

2.0 



124.1 ±Z? 

116.8 ±|;i 



1.69 
1.70 



2.80 

2.59 



— O 1U 
I 0.12 
±0.11 



2.86 ±°;|| 
2.38 ± ;2l 



195 +0.010 
— +08 

124 +" ■ our ' 
±0.005 

0.098 ±8 £11 



212.9 ±2 



9.0 
10.1 
(i.O 




262.9 ±lfl 
279.8 ±}g;f 



7 cq +0.81 
±0.76 
801 ±8:30 



5.44 
6.55 



1.6S 
1.11 

0.77 
0.08 



_o:ooo 

c 0.011 



0.140 l" "" 



Vikh 
Poly 



146.4 ±|;| 
138.3 ±H 



2.94 
2.72 



0.21 
25 
0.17 
0.17 



2.77 ±°;|i 0.106 ±8;<j" 
, 00 +0.30 n 117 4-0-008 



285.1 ±115 
278.6 ±}? o 



6.80 ±H 



4.09 ±S g2 0.170 ±S-S3I 



3.81 ± n 



o-i78 ±g; 024 



Vikh 
Poly 



213.3 ± 
189.1 ±g|f 



11.2 
7.2 
5.1 



3.31 
2.84 



3-59 ±°;| 
2.50 ±Q 20 



092 +(1005 
U.U»<! ±1) 008 

1) 1 1 8 -"-oo« 

0.113 ±0.005 



442.5 ±%% 
427.7 ±11J 



6.41 ±1 
5.79 ±9 



0.137 ±X 



Vikh 
Poly 



143.9 ±?;J 
138.4 ±f§ 



3.02 

2.87 



{823 
0.21 
O. 18 
0.13 



0.109 ±g;8ig 
0.116 ±§;i!g? 



_21.0 
-18.1 
L 13.6 
t 14.1 



4.26 ± 
4.65 ± 



1.02 
0.70 
0.88 
0.G4 



Vikh 
Poly 



Vikh 
Poly 



132.8 ±§;jj 

120.5 ±1;| 



+ 0.37 
±0.80 
+ 0.15 
±0. 15 



1.07 ±\\ 



2.75 ±S 



0.143 ±r ( 'Ai 

121 + 0-010 
u.±z± =t= n .009 



-n.014 

I] 1 50 -»3JO!) 
0.15O ±0.008 



211.0 ±i 



273.6 ±? n S s 



4.38 ± n l 



6.69 ±; : p 



o.i65 ±H:^| 

163 +0 024 
U.1.03 ±0.018 

0.164 ±0.084 
0.168 ±0;8;« 



Vikh 
Poly 



Vikh 
Poly 



Vikh 
Poly 



91.9 ± 
74.6 ± 



10.3 
13.5 
6.0 
5.4 



0.61 
0.36 



2.72 
2.30 



0.18 
0.28 
818 
0.00 



0.69 ±8;1| 



0.68 ±°;gjj 
0.66 ±» : 3fl 



3.31 ± 
2.34 ± 



1.41 

1.08 
8.88 
0.68 



0.089 ±g;8g| 
0.099 ±»;«S 



0.069 ±0;01| 



0-082 ±o : o;i 
0.098 ±0;0;i 



273.7 ±S;S 
319.0 ±tH 



4.14 ± 
5.56 ±; 



8.30 
8 40 
'.46 



E 0.45 
1 85 +030 
2.26 ± 



7.10 ± 
7.80 ± 



0.28 
0.66 

o'.Ya 



6.39 ±1;1| 



0.215 ±0;021 

0.190 ±s n ;g 



A 1 ^1 +088 8. 
0.131 ±0.025 
n 1 OS +0.029 
0.128 ±0.025 

0.112 ±0;02J 
0.078 ± -» 
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Appendices 



A APPENDIX: MODEL FITS TO 
TEMPERATURE PROFILES 

Temperature profiles for all the cluster s with the 
IVikhlinin et all |2006t ) model (blue) and the iBulbul et all 
model (red) are shown in Figure [121 Table [131 shows 
the radius out to which the temperature was measured and 
the estimated r^oo value for each cluster. Out of the 7 clus- 
ters, 5 clusters have temperature profiles which extended to 
^ 85% of rsoo- The temperature profiles of Abell 1758 and 
Abell 611 only reached 73% and 63% of rsoo, respectively. 
Both of these clusters suffered from solar flares, which did 
not allow the determination of the temperature profile out 
to rsoo. The effects of extrapolating the temperature profiles 
for these 7 clusters are not believed to significantly affect 
mass measurements. 
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Table 13. Extrapolation of Temperature Profiles out to rsoo 



Cluster 


Max Radius of T(r) 
(arcsec) 


»"500 
(arcsec) 


Per cent of rsoo 

(96) 


A1758 


240 


330 


73 


A1914 


420 


440 
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Figure 12. Temper ature profiles for all clusters using the 
IVikhlinin et al.1 (2006) model (blue) and the lBulbul et al.l [BoidD 
model (red). The solid lines show the best-fit values, and the 
hatched region is the 68.3% confidence interval. 
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